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Abstract: We investigate the cosmological properties of an 'induced gravity' brane sce- 
nario in the absence of mirror symmetry with respect to the brane. We find that brane 
evolution can proceed along one of four distinct branches. By contrast, when mirror sym- 
metry is imposed, only two branches exist, one of which represents the self-accelerating 
brane, while the other is the so-called normal branch. This model incorporates many of 
the well-known possibilities of brane cosmology including phantom acceleration (w < —1), 
self-acceleration, transient acceleration, quiescent singularities, and cosmic mimicry. Sig- 
nificantly, the absence of mirror symmetry also provides an interesting way of inducing a 
sufficiently small cosmological constant on the brane. A small (positive) A-term in this 
case is induced by a small asymmetry in the values of bulk fundamental constants on the 
two sides of the brane. 
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1. Introduction 

Extra-dimensional 'braneworld' models have attracted considerable attention in recent 
years. This is partly due to the fact that superstring/M-theory can only be consistently 
formulated in a universe which has more than four space-time dimensions. A distinctive 
feature of this theory is that it allows some of the extra dimensions to be large and even 
infinite thereby accommodating the braneworld scenario [1]. It is now well known that 
braneworld cosmologies can display quite distinctive behaviour which departs from that in 
general relativity either during early or late times [2] . The former can modify standard in- 
flationary predictions for primordial fluctuations while the latter can cause late-time accel- 
eration. This latter class of models will form the focus of our present paper. An abundance 
of recent cosmological observations points to a universe which is currently accelerating [3]. 
Although the source of cosmic acceleration remains unknown, most observations are well 
described by a small cosmological constant A/8irG ~ 1CT 47 GeV 4 . Such a small value for 
the A-term is difficult to explain within the context of standard field theory, which typically 
predicts a value for A which is several orders of magnitude larger than what is indicated 
by observations [4]. 

In this setting it is natural to ask whether cosmic acceleration could arise via an 
infrared modification of gravity at large distances. A famous example of such a scenario 
is the Dvali-Gabadadze-Porrati (DGP) model [5], which has been extensively discussed in 
many papers [2]. The present work shall examine a family of braneworld models which 
contain the DGP scenario as a subclass. As we shall show, this family has many interesting 
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features including that of phantom acceleration (w e g < — 1) and the possibility of inducing 
a small cosmological constant on the brane through bulk effects. 

Our braneworld scenario will contain a single large extra dimension. Two possibilities 
of principle exist in this case: either the bulk space is constrained to be symmetric with 
respect to the Z2 group of reflections relative to the brane, or such a symmetry is not 
imposed. The case where the bulk is symmetric is equivalent to the geometrical setting in 
which the brane is just a boundary of the bulk space; this is one way of justifying this mirror 
symmetry. An embedded brane without the Z2 symmetry is, however, a more general case 
with rich possibilities for cosmology. In spite of a considerable number of papers on this 
class of braneworld models (see, e.g., [6, 7, 8, 9]), some of these possibilities were either 
not noted or insufficiently studied. They form the subject of the present work. 

We consider a braneworld model described by the following simple yet generic action, 
which includes gravitational and cosmological constants in the bulk (B) and on the brane: 



Here, TZi is the scalar curvature of the five-dimensional metric g l ah on Bi, i = 1,2, the two 
bulk spaces on either side of the brane, and R is the scalar curvature of the induced metric 
hgj, on the brane. The quantity Ki = K l ab h ab is the trace of the symmetric tensor of extrinsic 
curvature K l ab of the brane in the space Bi. The symbol L(h a b, <t>) denotes the Lagrangian 
density of the four-dimensional matter fields <\> the dynamics of which is restricted to the 
brane so that they interact only with the induced metric h a b- All integrations over Bi and 
over the brane are taken with the corresponding natural volume elements. The symbols 
Mi, i = 1,2, and m denote the Planck masses of the corresponding spaces, Aj, i = 1,2, 
are the five-dimensional cosmological constants on the two side of the brane, and a is 
the brane tension. We shall focus on the asymmetric case with Ai 7= A2 and M\ 7^ M2 
which appears to be preferable from a string-theory perspective. For instance, the dilaton 
stabilised in different vacuum states on adjacent sides of the brane would lead to an effective 
five-dimensional theory with Mi 7= M2. The string landscape is likely to favour Ai 7= A2, 
which also occurs in domain wall scenarios. (The Randall-Sundrum (RS) [10] and Dvali- 
Gabadadze-Porrati (DGP) [5] models can be derived from (1.1) when mirror symmetry is 
respected.) 

In the absence of Z2 symmetry, cosmological evolution of the brane is described by 



where p is the total energy density of matter on the brane, Q = ±1, i = 1,2, correspond 
to the two possible ways of bounding each of the bulk spaces Bi, i = 1,2, by the brane. 
We classify the resulting four branches according to the signs of Ci and C2 as (++), (H — ), 



1 Note that, in the case of Z2 symmetry, there are only two ways of bounding the bulk by the brane, and 
these were called brane 1 and brane 2 in [11]. Of these, brane 2 contains the self-accelerating DGP brane as 
a subclass, while brane 1 can lead to phantom acceleration. 





(1.2) 



(-+), or {—)} 
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In the limit of Zi symmetry, the branches ( ) and (++) become the normal branch 

(brane 1) and self-accelerating branch (brane 2), respectively. The other two so-called mixed 
branches are characterised by C1C2 = — 1. As we shall show in this paper, cosmology on 
these branches can be quite novel. For instance, a small asymmetry in the values of the 
bulk constants on adjacent sides of the brane can induce a cosmological constant on the 
brane itself. On the other hand, in the limit when the bulk constants become exactly equal, 
one obtains stealth branes [7] that obey the general-relativistic internal equations but do 
not affect the metric of the bulk space. 

In this paper, we study the implications of (1.2) for a spatially flat universe (k = 0) 
without dark radiation (Cj = 0, i = 1,2). Equation (1.2) then simplifies to 



3m 2 m 2 ^— ' 1 1 V 6 3m 2 ^ £j 

i=l,2 i=l,2 

where we have introduced the fundamental lengths 



m 2 I 6 



£ * = wr A * = v-v i=1 ' 2, (L4) 

assuming negative values of the bulk cosmological constants. 

Note that (1.3) can be rewritten in terms of an effective cosmological constant, A e g, 

as 

where 

A e ff cr , (j / 2 2 



1=1,2 

which will be useful to us when we study the cosmological properties of this braneworld 
later in this paper. A pictorial representation of the branches described by (1.3) is given 
in the appendix. 

In this paper, we consider the late-time evolution of the universe, in which the energy 
density p is dominated by matter with the equation of state p = 0. Then, introducing the 
cosmological parameters 

° m = 3^f' "* = 3^H*' ^=^ 2 > ^ = A- 2 # - 2 , (1.7) 

where po and #0 are the current values of the matter density and Hubble parameter, 
respectively, we rewrite (1.3) in terms of the cosmological redshift z : 



H 



2/ 



h2 ( z ) = -T72 1 = fi m(i + zf + n a + V Q^hJh?{z) + n Xi . (1.8) 



^0 i=1>2 

This equation implicitly determines the function h(z), and explicitly the inverse function 
z(h). Note that the dimensionless cosmological parameters are related through the con- 
straint equation 

n m + n a +J2 dvWiV 1 + to)* = 1 . (1.9) 

i=l,2 

We now proceed to describe some features of braneworld cosmology without Z2 sym- 
metry. 
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2. Induced cosmological constant on the brane 



One way of accounting for cosmic acceleration within the framework of braneworld theory 
with mirror symmetry was suggested by Dvali, Gabadadze and Porrati [5]. An extension 
of this model to the case when mirror symmetry is absent is obtained by setting to zero 
the cosmological constants on the brane and in the bulk, so that a = 0, Aj = 0, i = 1, 2. 
The expansion law (1.3) then simplifies to 

i=l,2 1 



which evolves to a De Sitter limit at late times 



lim fr(*) = flD S = j.> (2 ' 2) 



1=1,2 ^ 

provided Yl Ci/^i is positive, which is true for branches (++) and (H — ), provided Mi > 

i=l,2 

M2 in the latter case. If m is of the order of the planck mass Mp ~ 10 19 GeV, then the values 
of Mi ~ 100 MeV can explain the observed cosmic acceleration. [The self-accelerating DGP 
model corresponds to the (++) branch with t\ = £2-] 

The absence of mirror symmetry provides a new avenue for this mechanism. Specifi- 
cally, the observed cosmic acceleration can be produced on one of the mixed branches with 
arbitrarily high values of the bulk Planck masses Mi and M2, provided these values arc 
sufficiently close to each other. If < AM = Mi — M2 <C Mi, then we have 

M, 3 - Mf 3M, 2 . . , 

H DS = 1 2 2 ~ -tAM 2.3 

on the (H — ) branch, and, by adjusting the value of AM, one can always achieve an 
observationally suitable value of Hds ■ For example, if Mi , M2 ~ m, then one needs AM ~ 
H . 

The previous model gave one example of late-time acceleration in the absence of the 
(brane) cosmological constant. We now derive another model with the same property but 
with a more flexible assumption A« / 0. Setting a = and A« / in (1.3) leads to 



i=l,2 

which evolves to a different De Sitter limit, expressed by the equation 



i=l,2 1 

where the length scales ii and Aj are defined in (1.4). 

It is interesting that a tiny asymmetry between the two bulk spaces can lead to a small 
cosmological constant being induced on the brane. Provided the bulk parameters Mi and 
M2 as well as Ai and A2 are close to each other, a neat cancelation on the right-hand 
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side of (2.5), which occurs for £1(^2 = — 1, leads to a small value of i^DS- Remarkably, this 
can happen even for very large values of the bulk constants. In particular, assuming that 
Aj <C -f^DS' we nave 

1 1 



4Ai e 2 x 2 (2,6) 

for one of the mixed branches. Thus, for bulk parameters of the order of a TeV, Mj ~ 1 TeV, 
Aj ~ 1 TeV^ 1 , we recover the current value of the Hubble parameter (i^os ~ Hq) provided 

Mi - ^ 2 A 2 | 1/2 - 1(T 13 TeV^ 1 - 1(T 30 cm . (2.7) 

Equations such as (2.3) or (2.6), (2.7) certainly represent fine tuning, with a tiny 
difference between bulk parameters only slightly breaking the smoothness of the metric 
across the brane. 2 In the limit of exact equality of the bulk constants on the two sides of the 
brane, the branches with (1(2 = — 1 describe a smooth bulk space, and the brane approaches 
the limit of a stealth brane [7] , evolving according to the usual Einstein equations without 
affecting the bulk space. 

As we shall see in the next section, the cosmological scenario with an induced cosmo- 
logical constant is distinguished by a property called cosmic mimicry which has interesting 
observational signatures. 

3. Braneworld expansion can mimic ACDM 

For large values of the bulk parameters, we encounter the phenomenon of cosmic mimicry 
which, in the context of Z2 symmetry, was described in [12]. Note that, during the radiation 
and matter-dominated epochs, the expansion of the universe follows the general-relativistic 
prescription 

where a/m? plays the role of the cosmological constant on the brane. However, at very 
late times, cosmic expansion gets modified due to extra-dimensional effects. Indeed, if 
Aj <C Hq 1 , then the square root in the last term of (1.3) can be expanded in the small 
parameter X 2 H 2 at late times, and the braneworld expands according to ACDM, namely 

* = ^ + I (3.2) 

with 

^(V , ,„-• I I - T ijfi ] CI.:!) 




*d^ESII'-E|! ■ (3.4) 



2 Perhaps, the small asymmetry in the fundamental constants characterising the bulk can be explained 
by the presence of the brane itself. For instance, the presence of the brane could lead to a small difference 
in the quantum contribution to the effective action of the bulk on its two sides, inducing slightly different 
bulk constants. 
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Note that both G and A are independent of time. Equations (3.2)— (3.4) have important 
ramifications. They inform us that the 'bare' value of the cosmological constant on the 
brane, a, is 'screened' at late times by extra-dimensional effects resulting in its effective 
value A. Thus, the early-time and late-time values of the cosmological constant are likely 
to be different, and this makes our model open to verification. 

Also note that one can have A / even if a = 0. Then, a small A-term can be induced 
during late-time evolution on the brane solely by extra- dimensional effects, as pointed out 
in the previous section. The mechanism by which the induced A-term becomes relatively 
small consists in a compensation of two potentially large terms with opposite signs in 
equation (3.4). Specifically, for small values of \ (which correspond to large values of Aj) 
such that Xi/£i <C 1, in the case a = 0, we have, approximately, 



which is another form of the result (2.6) for one of the mixed branches. What is remark- 
able here is that a positive cosmological constant on the brane can be sourced by bulk 
cosmological constants which are negative. 

From (3.2), (3.3) we also find that the effective gravitational constants during the early 
and late epochs are related by a multiplicative factor 



which can be larger as well as smaller than unity, depending on the braneworld branch. 
This factor will be closer to unity for the mixed branches (C1C2 = _ 1) than it is for the 
usual branches (C1C2 = 1) which survive in the case of Z2 symmetry. 

Focussing on the important case where a = and the effective four-dimensional cos- 
mological constant is induced entirely by five-dimensional effects, we find that, at redshifts 
significantly below the mimicry redshift 




(3.5) 




(3.6) 





(3.7) 



the brane expansion mimics ACDM 



h 2 (z) = n m (i + z f + n A , 



z <C Zj 



(3.8) 



with 'screened' values of the cosmological parameters: 




(3.9) 




-1 



(3.10) 



-6- 



On the other hand, from (3.2) it follows that, at high redshifts, the universe expands as 
SCDM 

h 2 {z) = Q m (l + z) 3 , z^z m . (3.11) 

An important distinguishing feature of this model is that the (screened) matter density, 
f2 m , inferred via geometrical tests based on standard candles and rulers, may not match 
its (bare) dynamical value Cl m . This allows cosmic mimicry to be distinguished from other 
cosmological scenarios by means of the Om diagnostic suggested in [13]. The fact that brane 
expansion also follows different laws at low and high redshift provides another important 
observational test of this model. 



4. Phantom branes 



In the presence of Z2 symmetry, the brane 1 branch of the generic model (1.2) exhibits 
phantom- like behaviour [11] which is in excellent agreement with observations [14] (see 
also [15]). Let us see whether this behaviour persists when mirror symmetry is absent. 
Note first that the condition for phantom acceleration, w(z) < —1, where 



w(z) 



2q(z) 



1 



3 [1 - n m ( z )\ ' 

has two equivalent formulations: 



?(*) 



dlogH(z) 
dlog(l + z) 



n-W > l d,osH{z) 



3dlog(l + z) 



Sl m (z) 



and A c g > , 



h 2 (z) 



(4.1) 



(4.2) 



where A e g- is the effective cosmological constant in (1.6), and differentiation is carried out 

with respect to the physical time variable. In the case of the ( ) brane (Ci = C2 = — 1), 

one has 



Aeff 



a 

3m 2 



E 



H 2 + A- 



(4.3) 



=1,2 



and we find immediately that A e g- increases with time when the expansion rate, H, de- 
creases. It is also quite clear that one (and only one) of the mixed branches will necessarily 
have a negative value of the sum term in (1.6), again exhibiting phantom behaviour. 

It is straightforward to verify that A e g > and H < on the two branches exhibiting 
phantom behaviour. Differentiating (1.3) and (1.6), we find 



H = — 



2-E 



< 0, 



A cff = 3HH 



for £ 



1,2 £i^H 2 + \: 



> 



-2 



(4.4) 
(4.5) 



i,2£i^H 2 + X- 2 



< 0. 
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Note that phantom models [16] with constant equation of state, w < — 1, are marked by 
A e ff > and super- acceleration: H > at late times. 3 This is related to the fact that 
the dark-energy (phantom) density in such models increases, as the universe expands, 
according to 

/Ophantom OC a^ 1+W ^ , W < -1 , (4.7) 

which causes the Hubble parameter to grow at late times, eventually leading to a Big-Rip 
singularity at which H diverges. By contrast, although the behaviour of our braneworld is 
phantom-like (w c g < — 1), the universe never super-accelerates since H < always holds. 
Furthermore, since H decreases during expansion, a Big-Rip-type future singularity which 
plagues phantom cosmology is absent in the braneworld. From the definition 

and property H < 0, we find q > —1. In fact, the deceleration parameter in our model 
always remains larger than the De Sitter value of q = — 1, approaching it only in the limit 
of t — ► oo. 



5. Transient acceleration 

An important property of this class of braneworld models is that the current acceleration 
of the universe need not be eternal. In other words, for a specific relationship between 
the fundamental parameters in (1.3), the acceleration of the universe is a transient phe- 
nomenon, and the universe reverts back to matter-dominated expansion in the future. 
Within the context of mirror symmetry, this scenario was called disappearing dark energy 
and discussed in [11]. In the absence of mirror symmetry, it was studied in [8] under the 
name stealth- acceleration (which should not be confused with the 'stealth brane' of [7]). 

Transient acceleration implies the property H — > in the asymptotic future, which 
requires the following condition to be satisfied: 



3^+Eix = ° =► ^+EWfyn* = o. (5.i) 



=1,2 i=l,2 



On the ( ) and (++) branches, this condition is realised with the following respective 

values of the brane tension: 



a 

3m 2 



(j^ + J^) ^ = ±(n/JV^+V^V^) ■ (5-2) 



3 It is easy to show that, in phantom models, the turning point H = occurs at 

1 + *. = -^= (I_^|1+hY /3H , w<-l. (4.6) 



a(t*) \ fi 
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On the new mixed branches (H — ) and ( — h), the required brane tension is smaller by- 
absolute value: 



a 

3m 2 



1A1 C2A2 



y/^^Ai - \A^ 2 ^A 2 • (5.3) 



=> n a = ± 

Under constraint (5.1), the cosmological evolution equation (1.8) becomes 

i — 1 9 \ 




(5.4) 



Condition (5.1) is necessary but not sufficient to speak about transient acceleration on 
a particular branch. A distinguishing property of a transiently accelerating brane is that 
q(z) — > 0.5 in the remote past (10 5 3> z 3> 1) as well as in the remote future (z — > — 1), 
reflecting the fact that the universe is matter dominated in the past and in the future, 
while, during the current phase, the deceleration parameter is negative, qo < 0. This last 
condition is realised only if the cosmological expansion law H 2 (p) is convex upwards to a 
sufficiently high degree. Specifically, in view of the second expression in (4.1), the condition 
qo < can be presented in the form 

dp 3p 

Looking at figures 3 and 4 in the appendix, one can see that this property can be realised 
only on two of the four branches: on the (++) branch and on one of the mixed branches. 

The expression for the current value of the deceleration parameter can be calculated 
by using the formula 

qo = ^ = -1. (5.6) 



One should note that only four out of five f2 parameters are independent in this expression 
because of the normalization condition /i 2 (0) = 1 applied to the evolution equation (5.4). 
In the ^-symmetric case, there remain only two independent f2 parameters. It is clear then 
that transient acceleration can be realised more easily in the ^-asymmetric case. This is 
illustrated in figures 1 and 2, which show the corresponding behaviour of the deceleration 
parameter q(z). 

In a transiently accelerating universe, cosmic acceleration is sandwiched between two 
matter-dominated regimes. A transiently accelerating braneworld clearly does not possess 
the Big Rip of phantom cosmology, nor even the event horizon of De Sitter space ! An 
in-depth study of this class of models [8] has revealed the existence of regions in parameter 
space which are stable (ghost-free). 

We have demonstrated that it is possible to construct braneworld models with transient 
acceleration. What is less clear is whether such transiently accelerating branches will pass 
key cosmological tests based on observations of high-redshift type la supernovae, baryon 
acoustic oscillations, etc. This important issue is open for further study. 



-9- 







" m = 0-3 


0.5 






1 o 










-0.5 




Symmetric braneworld 
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Figure 1: The deceleration parameter versus redshift is plotted for the (++) branch in the case 
of Z2 symmetry. The model has the parameters £t\ x = fl\ 2 = 2. We present plots for two different 
values of the matter density parameter: ri m = 0.3 and Q ra = 0.2. The sets of other parameters are 
calculated to be = f^ 2 = 1.21, fl a = —3.11 and = ilg 2 = 1.58, fl a = —3.56, respectively. In 
this case, an accelerated regime is not realised, although deceleration is significantly slowed down 
at the present cosmological epoch. 
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Asymmetric braneworld 



1 2 3 4 5 



Figure 2: The deceleration parameter versus redshift is plotted for the (H — ) branch in the case 
of absence of Z2 symmetry. The model has the parameters — 2 and il\ 2 = 2.1. We present 
plots for two different sets of values of the remaining two independent parameters: (flm^e^ — 
(0.3, 10000), which results in (fi^,^) = (9954.68, 3.16), and {il m ,il tl ) = (0.2, 5000), which results 
in (r^ 2 ,f2 CT ) = (4840.15, 0.82). Both plots show acceleration at the present cosmological epoch, 
which generically becomes more prominent for lower values of O m . 

6. Quiescent singularities 

A new feature of brane cosmology is a possible presence of quiescent singularities at which 
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the density, pressure and expansion rate remain finite, while the deceleration parameter 
and the Kretchman invariant, Riki m R tklm , diverge [17]. The universe encounters such a 
singularity in the future if a point is reached during expansion where the derivative of H 2 
with respect to p goes to infinity or, equivalently, where the derivative of p with respect to 
H 2 vanishes. Using (1.3), we can express this condition as the existence of a positive root 
H 2 of the equation 

i=i )2 iJm + K 2 



and a quiescent singularity is approached as H — > H s . At this moment, expansion formally 
ceases, and one cannot extend the classical evolution of the brane beyond this point. Such 
a singular point obviously exists on the (++) branch if and only if 



and it is reachable on this branch if the brane tension a is sufficiently negative: 

i<^ 2 -E7^<0, (6.3) 



3m 2 s ^ li 
i=l,2 



or, equivalently, 




<0. (6.4) 



Condition (6.1) may or may not be realised on the mixed branches. For example, in 
the simplifying case £\ = £2 = £, condition (6.1) is realised on the mixed branch (+— ) 
provided Ai > A2. One can show that the values of the parameters £i, Aj, i = 1,2, in 
principle can be chosen so that equation (6.1) has positive roots for three branches (++), 
(H — ) and ( — h). To achieve this, one only needs to satisfy the conditions l\ > £2 and 
Ai^2 > A2^i and choose sufficiently small values of £\, £2- 

For a graphical presentation of the reasons for the existence of quiescent singularities, 
the reader can look into the appendix. As in the case of mirror symmetry, quantum effects 
may play an important role in the vicinity of a quiescent singularity [18]; see also [19]. 

We also note that, in the case of mirror symmetry, realisation of quiescent singularity 
requires either negative brane tension or positive bulk cosmological constant (both condi- 
tions are suspicious from the viewpoint of possible instabilities). However the quiescent 
singularity can easily be realised without these assumptions in the asymmetric case on a 
mixed branch. 

The presence of a quiescent singularity in the future of the cosmological evolution does 
not threaten the past cosmological scenario. Therefore, this issue, just like the issue of Big 
Rip of phantom cosmology, is mainly of academic interest. Here, we only wish to point 
out that the possibility of quiescent singularity can be realised rather easily in braneworld 
theory in certain domain of its parameters without any additional ingredients (such as 
phantom fields, which lead to Big Rip singularities). 
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7. Discussion 



In this paper, we derive the expansion laws for an induced-gravity brane model in the 
general case where mirror (Z2) symmetry of the bulk space with respect to the brane may 
or may not be present. We find that, depending upon the choice of brane embedding, 
cosmological expansion on the brane can proceed along four independent branches, two 
of which survive in the case of Z2 symmetry. An important property of this class of 
models is that the four-dimensional gravitational and cosmological constants are effective 
quantities derivable from five-dimensional physics. In this case, brane expansion mimics 
ACDM at low redshifts, but the 'screened' matter density parameter f2 m does not equal 
its bare (dynamical) value f2 m . This opens a new avenue for testing such models against 
observations (see [12, 13] in this respect). Another important property of these models 
would be the growth of density perturbations which is likely to differ from ACDM (see 
[20], [21] and references therein). This issue lies beyond the scope of the present paper but 
we may return to it in the future. Braneworld models can be phantom-like and also exhibit 
transient acceleration. Thus, brane phenomenology, with its basis in geometry, provides 
an interesting alternative to 'physical' dark energy scenario's such as quintessence. 

The stability issues of the class of braneworld models without Z2 mirror symmetry 
were studied in [8, 9]. It is notable that ghost-free settings of the model with transient ac- 
celeration (and phantom-acceleration) appear to exist [8]. On the other hand, the analysis 
of the recent paper [9] reveals the presence of ghosts on a background with a De Sitter 
vacuum brane on the three branches (++), (H — ), ( — h) (i.e., which have at least one '+', 
so that the bulk at least on one side of the brane has 'infinite volume' in terminology of 
[9]). Whether this situation is critical for the cosmology under investigation remains to 
be seen. In this connection, it should be noted that the (++) branch, surviving in the Z2 
symmetric case, contains a ghost and is, therefore, linearly unstable [22]. On the other 

hand, the ( ) branch is ghost-free in the Z2 symmetric case and, apparently, also in the 

general case without Z2 symmetry. As we have demonstrated in this paper, this branch is 
responsible for 'phantom acceleration' (wo < — 1); also see [11, 14, 15]. 
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Appendix 

A. Visual representation of brane evolution 

The variables {X,Y}: 

X^f^L-H 2 , Y = H 2 , Ptot = p + a, (A.l) 
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allow us to rewrite equation (1.3) in the form 



Ci 'Y + X7 2 



(A.2) 



i=l,2 



which has a convenient visual interpretation. Equation (A.2) describes four branches in 
the physically restricted range Y > 0, with the symmetry of reflection with respect to the 
Y axis. If there exists a positive root Y c of the right-hand side of (A.2) for = — 1, then 
the two mixed branches intersect each other at the point 



Yc. = Hi 



(A.3) 



The condition for the existence of this intersection point is that the constant on the right- 
hand side of (A.3) be positive. 

The four branches (with and without intersection) are shown in Figs. 3 and 4, respec- 
tively. 





Figure 3: Four branches described by Eq. (A.2) in the (X,Y) plane and in the (p to t,H 2 ) plane 
in the case where two of them intersect. The horizontal dotted line indicates the position of the 
H 2 = axis in the case A1A2 = 0. The region below this axis is nonphysical. 



The brane evolves along one of the four branches towards decreasing values of ptot 
(during expansion). Depending (in particular) on the value of the brane tension a, three 
distinct possibilities can arise: 

(i) The trajectory may reach the value of H = 0, after which the universe recollapses 
and evolves along the same branch in the opposite direction. This happens when the 
value of ptot a t this point is greater than its minimum value a. 

(ii) The trajectory may asymptotically tend to either de Sitter space or the Minkowski 
universe with the minimum value ptot = °~- The second possibility occurs when 
the minimum value of ptot = c is exactly the point where the corresponding graph 
crosses the axis H 2 = 0, which, therefore, requires some amount of fine tuning. This 
possibility can be realised as transient acceleration. 
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Figure 4: Four branches described by Eq. (A. 2) in the (X,Y) plane and in the (p tot ,H 2 ) plane in 
the case of absence of intersection. The horizontal dotted line indicates the position of the H 2 = 
axis in the case A1A2 = 0. The region below this axis is nonphysical. 

(iii) The trajectory may end in a quiescent singularity at a finite value of H . This happens 
when the critical minimum point of ptot on the evolution curve is reached, and if this 
value of pt t is greater than its minimum value a. The reasons for the existence of 
quiescent singularities can be seen from the right panels in Figs. 3 and 4. They occur 
at the points of infinite derivative dH 2 /dptot- 
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